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ABSTRACT 

It has been shown before that the high mass-to-hght ratios of ultra compact dwarf galaxies 
(UCDs) can be explained if their stellar initial mass function (IMF) was top-heavy, i.e. that 
the IMF was skewed towards high mass stars. In this case, neutron stars and black holes would 
provide unseen mass in the UCDs. In order to test this scenario with an independent method, 
we use data on which fraction of UCDs has a bright X-ray source. These X-ray sources are 
interpreted as low-mass X-ray binaries (LMXBs), i.e. binaries where a neutron star accretes 
matter from an evolving low-mass star. We find that LMXBs are indeed up to 10 times more 
frequent in UCDs than expected if the IMF was invariant. The top-heavy IMF required to 
account for this overabundance is the same as needed to explain the unusually high mass-to-light 
ratios of UCDs and a top-heavy IMF appears to be the only simultaneous explanation for both 
findings. Furthermore, we show that the high rate of type II supernovae (SNII) in the star-burst 
galaxy Arp 220 suggests a top-heavy IMF in that system. This finding is consistent with the 
notion that star-burst galaxies are sites where UCDs are likely to be formed and that the IMF of 
UCDs is top-heavy. It is estimated that the IMF becomes top-heavy whenever the star formation 
rate per volume surpasses 0.1 Mq yr~^ pc~'^ in pc-scale regions. 



Subject headings: galaxies: dwarf — galaxies: star burst 
function, luminosity function — stars: neutron 



galaxies: stellar content — stars: mass 



1. Introduction 

The stellar initial mass function (IMF) quanti- 
fies the distribution of stellar masses in a newly 
born stellar population. Together with the de- 
pendency of stellar evolution on stellar mass and 
metallicity, as well as the rate at which stars are 
formed in the Universe, the shape of the IMF de- 
termines the chemical evolution of the Universe 
and how its stellar content changes with time. The 
shape of the IMF also has important implications 
for the evolution of star clusters. Thus, knowing 
the shape of the IMF is crucial for a broad variety 
of astrophysical problems. 

Resolved stellar populations in the Milky Way 
and its satellites support the notion that the IMF 



does not depend on the conditions under which 
star formation takes place, but that the stellar 
masses are distributed accor ding to a single IMF 
known as the canonical IMF ('Kroupa"200ll 120021 : 
Kumar et al.i .2008: Bastian_ct al. 201$"^ Ultra 



compact dwarf galaxies (UCDs) on the other hand 
provide evidence for the opposite notion, namely 
that the IMF varies and is top-heavy. 

These UCDs are stellar systems that have 
first bee n discovered in the Fornax galaxy cluster 
(jHilker et al. 1999). They have T^-band luminosi- 
ties between 10^ and some 10^ L(T), but half-light 



radii of only about 50 pc or less (jDrinkwater et al 



20031 : iMieske et ai1l2008f l. The confirmed UCDs 



are at distances where they cannot be resolved into 
stars with current telescopes, but constrains on 
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their stellar populations can be set by quantities 
derived from their integrated spectra. One such 
quantity are the dynamical mass-to-light {M/L) 
ratios of UCDs, i.e mass estimates based on the 



DM is predicted to gather on rath er large scales 
while UCDs are v ery compact (jGilmore et al 
,20071; iMurravl liooj ). However, in order to ex- 



density profile and the internal velocity dispersion 


of the UCDs (Hagegan et al. 


,20.05; Hilkcr et al] 


[ioOTt lEvstieneeva et all2007[ 


Mieske et al.ll200^ 



For a clear majority of the UCDs, the M/L ra- 
tios derived from their dynamics are higher than 
it would be expected if they were pure stellar 
copulations that formed with the canonical IMF 



populations tna t lormea witn tne canonical iM-b 
iHascgan et al. _ 2005t Dabringhausen et al. I2OO8I: 
iMieske et al.l I2OO8I ) . This has been taken as ev- 



idence for an IMF skewed tow ards high-mass 
stars ( Dabringhausen et al. 1 120091). i.e. a top- 
heavy IMF. The elevated M/L ratios of UCDs 
would then be explained by a large population 
of neutron stars and black holes (hereafter called 
dark remnants), because the age of the UCDs 
([Evstigneeva et al.ll2007t IChilingarian et al.ll2008l ) 



implies that all massive stars in them have com- 
pleted their evolution. 

It is plausible that the IMF in UCDs is skewed 
towards high-mass stars. Molecular clouds mas- 
sive enough to be the progenitors of UCDs become 
impenetrable for far-infrared radiation while they 
collapse and become a UCD-type star-cluster. In- 
ternal heating of the molecular cloud leads to a 
higher Jeans-mass in th em preferring the forma- 
tion of high- mass stars (jMurravl 120091 ). A molec- 
ular cloud can also be heated by an external 
flux of highly energetic cosmic rays originating 
from a local overabundance of type II supernovae 
incre asing the local Jeans-mass (jPapadopoulos 
2010[ ). With the young UCDs being very com- 
pact, also crowding of proto-stellar cores and 
their subsequent merging in young UCDs may 
lead t o an ov erabun dance o f high-mas s stars in 
them (jDabring hauscn et al Wcid ner et al 

201l[ ). 



However, the high M / L ratios of UCDs could in 
principle also be due to non- baryonic dark matter 
(D M), as was suggested by ICoerdt et al.l (|2008[ ) 
and lBau mgardt fc Mieskj ( 2008h . This is because 
dark remnants and non-baryonic DM would have 
the same effect on the M/L ratios of the UCDs, 
provided that a large enough amount of non- 
baryonic DM can gather within the UCDs. Note 
that that non-baryonic DM is an unlikely cause for 
the high M/L ratios of UCDs, since non-baryonic 



elude this possibility completely, the presence of 
a sufficient number of dark remnants has to be 
confirmed independently by a method that does 
not rely on the fact that dark remnants are non- 
luminous matter like non-baryonic DM. 

Such a method is searching for low-mass X-ray 
binaries (LMXBs) in UCDs. In these binary sys- 
tems, a dark remnant and an evolving low-mass 
star are orbiting around each other. The expand- 
ing outer atmosphere of the low-mass companion 
is accreted by the dark remnant. This matter pro- 
duces a characteristic X-ray signature. The num- 
ber of LMXBs depends on the number of NSs and 
stella r-mass blac k holes (BHs ) and thus on the 
IMF (|Verbunt fc H ut 1987; Verbimtll2003l ). This 
implies that stellar systems with a top-heavy IMF 
can be distinguished from stellar systems with the 
canonical IMF by an excess of LMXBs. 

The formulation of the IMF that is used 
throughout this paper is introduced in Sec. In 
Sec. ([3|), the LMXB-abundance in globular clus- 
ters (GCs) and UCDs in dependency of the IMF 
and this model is compared to observations. The 
type-II supernova rate in star-bursting galaxies 
in dependency of the top-heaviness of the IMF 
is discussed in Sec. (|H). It is found in Sec. ^ 
and Sec. (|H), respectively, that the UCDs and the 
star-bursting galaxy Arp 220 show indications for 
a top-heavy IMF. This suggests that the star for- 
mation rate per volume is perhaps the parameter 
that determines whether the IMF in that volume 
becomes top-heavy, as is argued in Sec. ([5]). Con- 
clusions are given in Sec. 

2. The initial stellar mass function 

A varying IMF can be formulated as 

i(m) =kk,m-'^\ (1) 

with 



ai = 1.3, O.K <0.5, 

- Mo 

m 

a2 = 2.3, 0.5 < — — < mtr, 

IVl© 

m 

as e M, mtr < ^TT- < rrimayi, 

Mq 
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where m is the initial stellar mass, the factors 
ki ensure that the IMF is continuous where the 
power changes and k is a normalization constant, 
^(m) equals if to < 0.1 or to > mmax, 
where TOmax is a function of t he star-cluster 
mass ( Weidner fc Kroupa 2006t Weidner et al. 
2OIOI ) and TOtr is the stellar mass at which 



the IMF begins to deviate from the canonical 
IMF. For TOtr — IM0, the formulation of the 
IMF used here is identi cal w ith the one used 

(|2009t ). so that results 



Dabringhausen et al 



found here for this choice of TOtr can be com pared 

(I2n09h . 



to results in iDabringhausen et al.l ()2009l ). For 
= a2 = 2.3 , Equa tion ([1]) is the canonical IMF 



(lKroupal[20Q]l l2002l) . For 03 < 2.3, the IMF is 
top-heavy, implying more intermediate-mass stars 
and in particular more high-mass stars. 

In the mass range of UCDs, TOmax is not set 
by the mass of the stellar system, but by the 
observed mass limit for stars, TOmax*- Thus, 
'TT-max = ™max* for all UCDs. The actual value 
of TOmax* is, however, rather uncertain: Estimates 



range from the canonical v a lue m 



(|Weidner fc Kroupal 120041: lOev fc Clark. 
to TOmax* ~ 300 Mp (ICrowthcr et al 



150 Mr, 



Baneriee et al.l 120111) . In this paper. 



2005) 



201C 



^ma 

but see _ 

Wmax* = 150 Mq is assumed, but note that as- 
suming TOmax* = 300 Mq instead would have little 
effect on the results reported here (see Section 13.31 
and Figure [7]). 

In the case of GCs and UCDs with LMXBs (see 
Section [3]), the observed luminosity, L, is known 
to originate from stars with masses to < 1 Mq. 
This is because their stellar populations a.re old 
(jEvstigneeva et alllioOTi IChilingarian et al.ll2008l ) 



and the more massive stars have already com- 
pleted their evolution. Being fixed by observa- 
tions, L should however not be changed when the 
IMF is varied. For the IMF given by Equation ([T]), 
this can be achieved by finding k from the condi- 
tion 

/ ^can{ni)m dm = IMq, (2) 

JQ.I Mq 

where ^can is the canonical IMF, i.e. = 2.3. 
With this normalization, the number density of 
stars with to < 1 Mq is the same for all values of 
as, since the normalization is set by the canonical 
IMF and is therefore not affected by variations of 
as- 

In the case of the SN-rate of Arp 220 (see Sec- 



tion m , the light used to estimate the star for- 
mation rate (SFR), i.e. the mass of the material 
converted into stars per time- unit, originates from 
stars over the whole range of stellar masses. With 
the SFR thereby given, we then normalize the IMF 
such that the SFR remains constant when the IMF 
is varied. For the IMF given by Equation ([T|) , this 
can be achieved by finding k from the condition 



^ {m)m dm — 1 M(; 



(3) 



0.1 M0 



With this normalization, the number density of 
stars with to < 1 Mq decreases with decreasing 
values of as, i.e. with increasing top- heaviness of 
the IMF. 

Stellar evolution and dynamical evolution turn 
the IMF of a star cluster into a (time-dependent) 
mass function of stars and stellar remnants; the 
star and stellar remnant mass function, SRMF. 
For a single-age stellar population, the connection 
between the IMF and the SRMF can be quantified 
by an initial-to-final mass relation for stars, TOrcm, 
which can be written as 



TOrc 



f TO 
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Mo' 
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M^ + «-^^4' 
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- Mo 
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(4) 

where TOto is the mass at which stars evolve 
away from the main seque nce at a given age 



(jPabringhausen et al.l 120091). UCDs typicall 



have ages of ~ 10 Gyr ([Evstigneeva et al 



icaiiy 
2007t 



Chilingarian et al.l 120081 ) . which implies TOto ~ 
1 Mo for them. In the present paper. Equation (|4]) 
is used to calculate how the mass of a modeled 
UCD depends on the variation of its IMF (see 
Section (3X3]). 

Note that Equation (|4]) reflects the evolution 
of single stars. In a binary system, the initial 
mass of a star that evolves into a black hole is 
expected to be higher, so that stars with masses 
up to 40 Mq may become NSs instead of BHs (cf. 
Brown et al.ll200l[ ). It is however of minor impor- 
tance in this paper whether a massive remnant is 
a NS or a BH. Both kinds of objects can become 



3 



bright X-ray sources by accreting matter from a 
companion star and BHs in such binary systems 
are actually detected b y excludi ng that they are 
NSs due to their mass (|CasaresH2Q07t ). Also the 
total mass of a GC or UCD is not strongly af- 
fected by the mass-limit between NSs and BHs. 
Using Equation (refeq:IMF) with mtr = 1 M0 and 
Equation (jH) with mto = 1 M©, the total mass of 
NSs and BHs is 4.2 per cent of the total mass of 
the stellar system for = 2.3 (canonical IMF) 
and 79.9 per cent for 03 = 1. These numbers are 
altered to 3.8 per cent of the total mass of the 
stellar system for 03 = 2.3 and 75.0 per cent for 
a3 = 1 if the transition from NSs to BHs is shifted 
from 25 Mq to 40 Mq. 

3. The LMXB-abundance in GCs and 
UCDs 

3.1. Some properties of GCs and UCDs 



Fo r a number of GCs and UCDs, data (|Mieske et al 
20081 ) on U-band luminosity (iy), dynamical 
mass (Mdyn) and effective half-light radius (rh) 
are available. These data suggest a transition at 
Lv ~ 10® Lq, since the rh and dynamical M/L 
ratios of objects with Ly < 10® Lq appear to be 
independent of Ly, in contrast to objects with 
Lv > 10® Lq (see Fi gures [1] and [2|) . This moti- 
vates to consider the objects with Ly < 10® L© 
as GCs and those with Ly > 10® L© as UCDs, 
even though stellar systems close to this transition 
could be assigned to either one of these classes 
iMicskc ct_al,, 2008j). 

Knowing rh and Mdyn of a stellar system allows 
to estirnate it s median two-body relaxation time 
(|Spitzeiill987l ). using 



^rh — 



0.234 



logio(A/dyn/Mo) 



G 



(5) 



the ACDM- model (see Figure |3]). Thus, the prop- 
erties of UCDs can be calculated from their IMF 
while considering the effects of stellar evolution, 
but without accounting for the effects of dynam- 
ical evolution. This means in particular that the 
SRMF of UCDs can be calculated from their IMF 
and Equation ([4]). Note that GCs, on the other 
hand, are subject to dynamical evolution, since 
their ages, tgc are also similar to th and thus 

2008[) on Lv and rh 



The data ( Mieske et al 



of individual GCs in the MW and in Centaurus 
A and UCDs in the Virgo-cluster are also useful 
for estimating an average rh, rh, as a function of 
Ly- GCs over the luminosity range from 10* L© 
to 10® L(7) do not show a luminosity-r adius trend 
(|McLaughlinll2000l: I Jordan et al.ll2005[) . The loga- 
rithmic average rh of GG is 



logio ( 



0.4314 



(6) 



([Jordan et al.l 120051 ) . Performing a lin ear least 



squares fit to data in lMieske et al. 1 (120081) on UCDs 
in the Virgo cluster leads to 



pc 



logio — = 1.076 log 



10 



Li 



10® L 



0.4314. 



wher e G is the gravitational constant (jDabringhausen et 
20081 ). The significance of tj-h lies in the fact that 
it sets the time-scale on which the structure of a 
self-bound stellar system is changed by the pro- 
cess of energy equipartition. If t > tj-h holds for a 
stellar system with r being its age, it can be con- 
sidered nearly unaffected by dynamical evolution 
and is thus only subject to stellar evolution. This 
is the case for UCDs, as i^h ?J th is valid for them, 
where th is the age of the Universe suggested by 



(7) 

Note that equality between Equations ([6]) and ^ 
at Lv = 10® L© was imposed as a secondary con- 
dition on the fit of Equation ^ to the data. This 
secondary condition reflects the fact that the rh of 
GCs are indistinguishable from those of UCDs at 
L w 10® L© (see Figures [1] and [J) . 

3.2. Modeling the LMXB-abundance in 
GCs and UCDs 

3.2.1. The origin of LMXBs in GCs and UCDs 

1 Tight binaries consisting of a dark remnant and 

-Sla, low-mass companion can have in principle two 
different origins: 

1. They can be primordial. In this case a tight 
binary of a high-mass star and a low-mass 
star have formed already in the star form- 
ing event. The high-mass star explodes in 
a supernova after a few million years leav- 
ing behind a dark remnant which can remain 
bound to its low-mass companion. 
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Fig. 1. — The effective half-light radii, rh, of GCs 
and UCDs. The circles show the sample of in- 
dividual GCs and U CDs from the compilation of 
Mieske et al. The dashed line is an es- 



timate of the average rh of GCs and UCDs (cf. 
Equations [5] and [T]) . The vertical dotted line sets 
the limit between objects that are considered as 
GCs and objects that are considered as UCDs. 
Note that the average rh indicated for GCs by 
the dashed line is lower than the average rh of 
the GCs shown in this figure. This is because the 
GCs shown here are mostly GCs of the Milky Way 
while the dashed line corresponds to the average 
rh of GCs in the Virgo-cluster. The GCs in the 
Virgo cluster tend to be more compact than those 
around the Milky Way. 



2. They have formed through encounters. GCs 
and UCDs are regions of enormously high 
stellar density ranging fr om 10 Mf?) p c~^ to 
10''' Mq pc~'^ (iDab ring hausen et al. 



l2008l) . 



Encounters between dark remnants and low- 
mass stars are therefore frequent and can 
lead to the formation of LMXBs due to 



20031) 



tidal capture (jVerbunt fc HutlllQSTt IVerbunt 



As these formation mechanisms are quite diff^er- 
ent it is expected that both processes would con- 
tribute differently to the LMXB content in GCs 
and UCDs. 
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Fig. 2. — The mass-to-light ratios (M/L ratios) of 
GCs and UCDs. The circles show the sample of 
individual GCs and UCDs from the compilation of 
Mieske et aTl (|200i). The vertical dotted line sets 
the limit between objects that are considered as 
GCs and objects that are considered as UCDs. 



There are however strong arguments against a 
significant contribution from primordial binaries 
to the LMXB content of GCs and UCDs: 

1. The number of LMXBs in GCs is strongly 
correlated with the encounter rate an d thus 

20051: 



clearly linked to it ([Jordan et al 



Sivakofl[ et al.ll2007D 



2. There are several hundred times more 
LMXBs per u nit mass in GCs than in the 
Galactic field (jVerbunt fc Hutlll987l ). The 
LMXBs in the Galactic field are LMXBs 
that probably evolved from primordial bi- 
naries, since they are in a low-density envi- 
ronment where encounters play no role and 
most probably formed in star clusters from 
which they were subsequently ejected. The 
strong excess of LMXBs in GCs therefore 
suggests that most JjMXBs in GCs form 
through encounters (iVerbunt & Hut 1987). 

The number of encounters relevant for the cre- 
ation of LMXBs, i.e. encounters where a NS can 



capture a low- mass star ([Verbunt fc Hutj llQSlT) , 
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tween NSs and red giant stars or interactions be- 
tween stars and existing binaries also play a role 
and can actually be even more important than 
tidal captures. However, this does not change the 
observational finding that th e number of LMXB s 
in GCs scales with T (e.g. Jordan et al. 2005[) . 
Therefore, F seems to be an adequate measure 
for the stellar dynamical processes that produce 
LMXBs in g e neral . It is moreover argued in 
Ivanova et al. that primordial binaries only 

make a small contribution to the total population 
of LMXBs in old GCs. 

There are thus strong observational and theo- 
retical motivations for the usage of F as a measure 
for how many LMXBs are expected in GCs and 
UCDs. 



Fig. 3. — The median two-body relaxation times, 
trh, of GCs and UCDs. The circles show the sam- 
ple of individual GCs and UC Ds from the com- 
pilation of .Mieske et al.l (|2008l ) . The dashed hor- 



izontal line indicates the age of the Universe, th, 
according to the ACDM-model. The vertical dot- 
ted line sets the limit between objects that are 
considered as GCs and objects that are consid- 
ered as UCDs. Note that i^h > m for UCDs. 
Thus, U CDs can be considered dy namically un- 
evolved ( Dabringhausen et al.ll20Q8 ) and they may 
therefore be considered as galaxies from a stellar 
dynamical point of view fjForbes fc Kroupa,201L) . 



can be written as 



Foe 



(8) 



where rins is the number density of NSs, rig is the 
number density of potential low-mass companion 
stars, Tc is the core radius of t he stellar syste m and 
a is the velocity dispersion (|Verbunti 120031 ). The 
potential companions to a NS in a bright LMXB 
are stars that come from a rather narrow mass 
range where stars of a given age leave the main- 
sequence. At this stage of their evolution, the stars 
expand rapidly, which makes a high accretion rate 
on the NS possible, which in turn leads to a high 
X-ray luminosity. 

A more recent study by llvanova et al. I (l2008l) 

revealed that tidal capture is not the only dynam- 
ical process relevant for the formation of LMXBs. 
Other dynamical processes like direct collisions be- 



3.2.2. The encounter rate in GCs and UCDs for 
an invariant SRMF 

If only a single, invariant mass function for stars 
and stellar remnants (SRMF) is considered for all 
stellar systems, then 



oc po 



(9) 



holds, where po is the central mass density. Equa- 
tion ([8]) can then be rewritten as 



F oc 



pI 4 



(10) 



by using Equation (0). 

In order to link the theory on LMXB-formation 
to the optical properties of observed stellar sys- 
tems, it is in the following assumed that the mass 
density of a stellar system follows its luminosity 
density. The structural parameters derived from 
the distribution of the light in the stellar system 
can then be translated directly into statements on 
the distribution of its mass, i.e. quantities that 
determine the dynamics of the stellar system. 

However, Tc is difficult to measure for GCs and 
UCDs at the distance of the Virgo cluster, as these 
stellar systems are barely resolved with current in- 
struments. The projected half-light radius Vh (and 
thus the half-mass radius under the assumption 
that mass follows light) is larger and therefore less 
difficult to retrieve from the data. For practical 
purposes, it is therefore useful to assume 



Tc oc Th 



(11) 



6 



and 



M 

Po oc -J, 



(12) 



where M is t he mass of the stehar system. The 
King profile (|KinsJll962l) with its three indepen- 
dent parameters (core radius, tidal radius and 
central density), is thereby simplified to a den- 
sity profile with only two independent parame- 
ters (half-mass radius and mass). The underly- 
ing assumptions are not necessarily true, and in- 
deed, not fulfilled for GCs in the Milky Way since 
[McLaughlin (200 0) finds that more luminous GCs 
tend to be more concentrated. However, regard- 
ing the conclusions on how the presence of bright 
LMXBs is connected to the optical properties of 
GCs in the Virgo cluster, these assumptions are 
unproblematic. Using t he same conce ntration for 
all GCs in their sample ISivakofF et al. 1 (2007) find 
that they ess e ntially come to the same results as 



Jordan et al.l ( 20041 ) . who use an individual esti- 



mate for the concentration of each GC in their 
sample. 

When dealing with UCDs, replacing Tc with rh 
is even advantageous. The time-scale on which the 
NSs gather at the centre of the UCD is given as 



^scg — 



■t. 



(13) 



where m « 0.5 Mq is the mean mass of stars, 
mns ~ 1-35 Mq is the mass of neutron stars and 
^cc is the core-collapse time of the UCD w ithout 



the NSs (ISpitzeijll987t iBaneriee et al.ll2010t ). If a 
Plummer sphere ( Plummeijll91ll ) is used as an ap- 
proximation for the densi ty profile of a stellar sys - 
tem, tec ~ 15 irh holds (jBaumgardt et al.ll2002l) . 
With ifh being of the order of a Hubble time for 
UCDs, Equation (IT3|) implies that the distribution 
of NSs in UCDs still follows the initial distribution 
of their progenitors. The volume relevant for the 
formation of LMXBs in a UCD is therefore bet- 
ter measured by than by Tc , provided its stellar 
population did not form mass-segregated. This is 
because rh represents the size of the whole UCD, 
whereas Tc represents the size of its centre. 

Thus, using Equations (jlip and ()12|) . Equa- 
tion pUj) can be transformed into 



Lh oc 



(14) 



If the stellar system is also in virial equilibrium, 

a oc Po' ''h oc — p-g- (15) 
holds. In this case. 



oc -33- (16) 



follows from Equations (|14p and (|15l) . In contrast 
to Equation ([5]), Equation has only two vari- 
ables (M, rh) instead of four (rig, rins, Af, f'c)- 

A further variable can be eliminated by re- 
placing individual values for rh by luminosity- 
dependent estimates for rh, such as Equations ^ 
and (O, and noting that the same SRMF for all 
stellar systems in question implies M oc Ly for 
them. This leads to 



Fh oc 



rl.5 



(17) 



or, more explicitly by using Equations ([5]), and ([7]), 
respectively. 



logio(rh) = l-51ogi 



106 L. 



for GCs (i.e. Ly < 10^ Lq), and 



logio(rh) = - 1.190 log 



10 



Li 



106 U 



+ A (18) 



+ A (19) 



for UCDs (i.e. Ly > 10^ L©). The constant A is 
the same in Equations (fT8|) and ([T9| . Note that 
the transition between Equations (fT8|) and (fT9|) is 
continuous due to the continuity of fh at Ly — 
106 Lq. 

3.2.3. Detecting a variable SRMF with LMXBs 

For investigating how Fh depends on the IMF, 
it is useful to consider the ratio between Fh as a 
function of and the Fh implied by some ref- 
erence IMF. This has the advantage that factors, 
which do not depend on the IMF, cancel. The ref- 
erence IMF is the canonical IMF in this paper; a 
choice that is motivated with the lack of dynam- 
ical evolution in UCDs (cf. Section |3TT|) . Using 
Equation (fT2|) thus leads to 



rh(a3) 



iM{a3 = 2.3) 



rh(a3 = 2.3) n„s(a3 = 2.3) 



M(a3) 



(20) 



7 



if it also assumed that the IMF varies only for 
stars with m > mto, so that also rig is constant. 
By this last assumption, the luminosity of the 
UCDs, which is given by observations, stays con- 
stant when the IMF of the UCDs is varied. The 
right side of Equation ([20)) can be calculated if the 
IMF is specified. In particular. 



"■ns(a3 = 2.3) J^mT ^can(m) dm ' 



(21) 



and 



M(a3 = 2.3) /oTmq "^'■cm ('^)€can (to) dm 



Jo.l Mq "^remv 



{rn)^{m) dm, ' 

(22) 

where the IMF is normalised using Equation ©, 
^can is the canonical IMF and mremijn) is given 
by Equation (jj]). Thus, Equation ([20)) quantifies 
how Fh changes in a stellar system (normalized 
with Fh for the canonical IMF) if the number of 
dark remnants and therefore the mass of the stellar 
system are changed, while its characteristic radius 
and the number of stars are kept constant. 

A difficulty is that Fh of a stellar system can- 
not be measured directly. However, the actual F^ 
of a GC or a UCD scales with the rate at which 
LMXBs are created (see Sec. 13. 2. 1")) , which is pro- 
portional to the probability P to form an LMXB 
above a certain brightness limit in a given time. If 
a sample of GCs or UCDs in a certain luminosity 
interval is given, a useful estimator for the aver- 
age P of these GCs or UCDs is the fraction /lmxb 
of them that have an LMXB above the brightness 
limit defined by the sensitivity of a given set of 
observations. Thus, 



/lmxb oc P cx F^, 



(23) 



where the exponent 7 accounts for the claims 
that the LMXB-frequency in GCs and UCDs may 
not be directly proportional to F or Fh, but to 
some power of F o r Fh (cf. Ijordan et al. 2004 : 
Sivakoffet ai1l2007h . 

If the SRMF of UCDs is indeed independent of 
luminosity, the /lmxb of UCDs in different Ly in- 
tervals should all roughly coincide with the predic- 
tion from Equation ()19p for an appropriate choice 
of the constant A. If however the /lmxb of at 
least one Ly interval is inconsistent with Equa- 
tion ([19)) for any choice of A, then this would be 



evidence for the SRMF changing with the lumi- 
nosity of the UCDs. This would imply that the 
IMF of the UCDs changes with luminosity, since 
the SRMF of UCDs is solely determined by stellar 
evolution, i.e. a process that does not depend on 
the size of the system (see Section [3TT)) . Note that 
the actual value of A in Equations ([TO)) and ([TO)) 
has no implications for the physical properties of 
the observed stellar systems: For a given sample 
of GCs and UCDs, A depends on the detection 
limit for an X-ray source or an arbitrarily chosen 
brightness limit above the detection limit. 

3.2.4- Data on the LMXB-frequency in GCs and 
UCDs 

In order to search for a dependency of the 
IMF in UCDs on their luminosity, we use data 
published in the up per left panel of figure (6) in 
Sivakoff et al.l (|2007l) . These data provide the frac- 
tion of globular clusters and UCDs, /lmxb, host- 
ing an LMXB in a given total z-band magnitude 
interval. 



The results of ISivakoff et al.l (|2007|) were ob- 
tained by combining two sets of data. 

First, HST images of 11 elliptical galaxies 
in the Virgo Cluste r were used, see Table 1 in 
Sivakoff etall (|2003). Ten of them are the bright- 
est galaxies observed ii i the course of t he ACS 
Virgo Cluster Survey (|C6te et al.l |2004| ). The 



eleventh one (NGC 4697) is a similarly bright 
galaxy that was observed bv lSivakofreTal1(l200i1) 
with nearly the same observational setup as in the 
ACS Virgo Cluster Survey. Using the obtained 
images, a large number of accompanying GCs and 
UCDs was identified around each of these galaxies. 
Second, ISivakoff et al.l (|2007l ) used archival 



Chandra Observatory X-ray observations of the 
same galaxies. The setup for the X-ray obser- 
vations var ied widely from gala xy to galaxy, see 
Table 2 in ISivakoff et al.l (|2007t) . which could in 
principle be problematic. 



Sivakoff et al. 



([2007D find however that the 
global properties of GCs and UCDs which contain 
a LMXB are largely unaffected by the varying de- 
tection limits for X-ray sources. Also note that 
the LMXB-frequencies in GCs are well explained 
by the encounter rates in them (see Section 13.31) , 
despite the different detection limits for X-ray 
sources. This suggests that the encounter rate 
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is indeed a good measure for the rate at which 
LMXBs of any X-ray luminosity are created. We 
therefore assume that a large number of GCs and 
UCDs with an X-ray source is indeed an indicator 
for a large number of dark remnants in them. 

The size of the z -band magnitude intervals in 
Sivakoff et al.l (j2007l ) is chosen such that each of 
them contains 27 GCs or UCDs with a detected 
LMXB. This corresponds to a total of at least 
100 GCs or UCDs in each of these intervals, since 
/lmxb < 0.2 in ah of them. Thus, /lmxb can 
be taken as a reliable estimator for the average P 
to form an LMXB in a GC or a UCD in a given 
z-band magnitude interval. 

For comparin g the data on the LM XBs in GCs 
and UCDs from ISivakoff et al.l (|2007l ) to the pre- 
diction for the LMXB-frequency in GCs and UCDs 
formulated in Equations ([TSl) and ([T9l) . z-band 
magnitudes have to be converted into Ly- For 
this purpose, z-band luminosities are calculated 
from z-band magnitudes with 



(24) 



where is the absolute z-band magnitude and 
Lz is z-b and luminosity in S olar units (cf. Equa- 
tion 1 in lSivakoff et al.|[2007h . Now note that the 
z-band M/L ratio of GCs in the Virgo-cluster ar e 



all close to ?s 1.5Mq/Lq^, (jSivakoff et al.ll200iD . 



which is essentially identical to the average V- 
band M/L ratio of the GCs in the Milky Way in 
Solar units (|McLaughlinll2000[ ). This implies that 
z-band and T^-band luminosities of GCs are ap- 
proximately identical in Solar units. We therefore 
assume Ly/LQy = L^/Lq^z in this paper. 
The data from figure (6) 



Sivakoff et al 



POOTI ) is shown in FigurelHwith the z-ba nd magni- 
tude intervals from Sivakoff et al.l ( 2007 ) converted 
into Lv intervals. Three of these intervals are at 
luminosities Ly > 10® Lq, so that the objects 
in them are UCDs (cf. Section [3TT|) . As the size 
of the intervals is chosen such that each of them 
contains 27 objects with an LMXB, 81 UCDs with 
an LMXB are considered here. The total number 



of UCDs in the sample from lSivakoff et al.l (|2007|) 
is about 400, as can be calculated from /lmxb in 
the according intervals. 

For practical purposes, it is useful not to dis- 
cuss individual values for the /lmxb of UCDs, but 
to replace them by a continuous function P(Ly). 



This function is obtained by performing a least- 
squares fit of a linear function to the values for 
/lmxb in the Ly intervals with the UCDs, lead- 
ing to 

logi„(P) = alogio(Ly) + 6, (25) 

where the best fitting parameters a and b are given 
in Tab. P can be interpreted as an estimate 
for the average probability for UCDs with a given 
Lv to host a LMXB brighter than the detection 
limit. For a meaningful comparison between P 
and Fh at different values for Ly, A needs to be 
gauged. This is done by imposing that P{Lv) = 
V\i{Lv) for Lv = 10® Lq^. The motivation for 
choosing this condition to fix A is that the stellar 
populations of systems with this luminosity should 
be nearly unaffected by dynamical evolution (cf. 
Section l3TT|) . while their M/L-ratios suggest that 
their IMF is canonical, in c ontrast to even more 
lumin ous stellar systems (cf. iDabr inghausen et al 
20091) . 



If the rate at which LMXBs are produced in 
GCs and UCDs is proportional to some power 7 
of the encounter rate in them, leading to P(Lv) oc 
Tl (cf. Equations [23l and [25]) . Equation ([SQ]) can 
be transformed into 



".ns(a3) 



PjLy)- _ 

rh(iy) nns(a3 = 2.3)' 



'M(a3 = 2.3) 



M{a3) 



(26) 



The left side of Equation (1261) is then expressed in 
terms of observable properties of UCDs and the 
right side only depends on as a free parameter 
once mtr is given. Equation (|26p can therefore be 
used to estimate the dependency of the IMF of the 
UCDs as a function of their observed Ly. Since 
A is chosen such that P{Lv)/T{Lv) — 1 for stel- 
lar systems that are assumed to have formed with 
the canonical IMF, P{Lv)/T{Lv) > 1 implies a 
top-heavy IMF and P(Ly)/r(iy) < 1 implies a 
bottom-heavy IMF. 

3.3. Results 

In order to test for an LMXB-excess and thus 
a top-heav y IMF in UCDs from the observational 
data from Isivakoff et ahl (2007) we now compare 
with theoretically expected LMXB-frequencies. 

The dynamical formation of LMXBs depends 
on the density of both dark remnants and low- 
mass stars (Equation |8]) . In denser star clusters. 
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Fig. 4. — The observed LMXB-frequency of GGs and UCDs in comparison to expected frequencies if the 
IMF was canonical. Plotted are the observed frequencies (squares) of GCs {Ly < 10^ io,v) and UCDs 
{Lv > 10^ L0,y) with LMXBs, /lmxb, in the Virgo galaxy cluster as a function of the V-band luminosity, 
Lv- Each data point contains 27 objects showing t he LMXB signa l . The data points are identical with the 
data points in the upper left panel of figure (6) in ISivakoff et al.l (|2007l ). except for a rescaling of z-band 
magnitudes to y-band luminosities. The three brightest data are based on w400 UCDs, i.e wl35 UCDs 
per bin. The dashed line shows the theoretically expected LMXB-frequency for an invariant canonical IMF 
with index ~ 2.3 assuming /lmxb cjc Fh (left panel) and /lmxb oc Fj^ (right panel). In either case, the 
theoretically expected LMXB-frequency is significantly too low for UCDs, while for GCs the theoretically 
expected LMXB-frequency matches the observed LMXB-frequency. The solid line is a fit through the UCD- 
regime (above 10^ L0^y). From it is derived the variation with luminosity of the IMF index such that 
this new model, based on a variable IMF, accounts for the observed /lmxb for Lv > 10^ Lq^v. 



close encounters are more frequent and the for- 
mation of an LMXB is more likely. GCs have 
a common half-mass radius of a few parsec in- 
dependent of their luminosity and their stellar 
mas s is on average pro portional to their luminos- 
ity ([McLaughlin It therefore follows from 
Equation ([T5| that LMXBs should be hosted pre- 
dominantly in high-mass GCs if their SRMF does 
not depend on their stellar mass. The dashed 
line in Figure |4] shows the theoretically expected 
LMXB frequency for a constant IMF calculated 
with Equations (fTS]) and (fT9| with A chosen 
such that these equations reproduce the observed 
LMXB frequency at Ly = 10^ Lq. The theoreti- 
cal prediction then matches the observations in the 
GC regime (i.e. Ly < 10^ L©), in ag reement with 



earlier studies on LMXB s in GCs (jJordan et al 
20051: iPeacock et al.|[2oTol ). 

At Ly « 10^ Lq, the transition luminosity 



from GCs to UCDs, both kinds of stellar sys- 
tems have the same half-mass radius (see Sec- 
tion 13. ip . However, unlike GCs, UCDs show 
a luminosity-radius relation such that they be- 
come less dense with increasing lumi nosity (cf. 
figure 4 in iDabringhausen et al.l |2008() . Conse- 
quently, Equation ([T51) predicts that the capture 
rate of late-type stars by dark remnants and thus 
the expected LMXB frequency decreases rapidly 
with increasing Ly-band luminosity if the SRMF 
is constant. Note that a constant SRMF in UCDs 
implies a constant IMF in them due to the lack of 
dynamical evolution in UCDs (see Section [Q]) . 

In Figure 01 the prediction from Equation (ITO)) 
for the LMXB frequency in UCDs with a constant 
SRMF is shown by the dashed line in the accord- 
ing luminosity range, where A is chosen such that 
Equation reproduces the observed LMXB fre- 
quency at Ly = 10^ Lq. Two cases are consid- 
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ered, namely /lmxb oc Th and /lmxb oc T^^. The 
second case is closer to th e dependency b etween 
/lmxb and Fh reported by ISivakoff et al ] 12007). 
The agreement between the theoretical prediction 
and and the observed frequency of LMXBs is in 
either case good for GCs. However, /lmxb oc r°-® 
seems indeed a better fit to the data than /t.m xb oc 
Th. Note that iMaccarone k Peacock! (|201l( ) find 



r oc Ff; on average for GCs in the Milky Way, 
which essentially means that the typical ratio be- 
tween F and Fh depends for these GCs on their 
mass. This is probably a consequence of the more 
massive GCs in the MW bei ng more concentra. ted 
than the less massive ones (|McLaughlinll200nl: cf 
Section [3.2.2p . and likely to be the case for the 
GCs in the Virgo cluster as well. 

For UCDs however, the observed LMXB- 
frequency strongly deviates from the theoreti- 
cal prediction for a constant SRMF. It is ob- 
served that 25 ± 5 per cent of the UCDs with 
« 5 X 10*^ Lq^ have a bright LMXB, while 
Equation ([T9| suggests a LMXB frequency of 
about 2 per cent at this luminosity for /lmxb oc Fh 
and a LMXB frequency of about 3 per cent for 
/lmxb oc r^-^. Thus, the expected fraction of 
LMXBs hosting UCDs is up to w 10 times smaller 
than observed if all UCDs had the same IMF. 

This discrepancy between the data and the 
model with an invariant (canonical) IMF and 
the data is highly significant. This cannot be 
explained with more dark remnants remaining 
bound to UCDs due to higher escape velocities. 
This is because the escape velocity from massive 
GCs is much higher than the escape velocity from 
light GCs, since the characteristic radii of GCs do 
no change with mass, but the encounter rate is 
nevertheless sufRcient for quantifying which frac- 
tion of them has a bright LMXB. 

The situation is more complicated with the 
finding that redder GCs and UCDs have more 
LMXBs than the blue ones, while brighter ob- 
jects (i.e. the UCDs in particular ) tend to be 



redder than the less luminous ones ( Mieskeetal 



20061 ). Taking color as an indicator for metallic 
ity leads to the interpretation that the LMXB 
frequency in GCs and UCDs does no t only depend 



on F or Fh but also on m etallicity pordan et al 



20041: ISivakoff et al.ll2007l) . Note that an increase 
of metallicity with luminosity and therefore mass 
of GCs is consistent with theoretical modeling. 



according to which more massive star clusters 
retain more processed (i.e. metal-enriched) gas 
which is turned into subseq uent stellar populations 
(|Tenorio-Tagle et al.ll20ol . 

Using metallicity (i.e. color) as a second pa- 
rameter besides Fh indeed allows a more precise 
modeling of the probability to find a LMXB in a 
given GC or UCD than when Fh is assumed to be 
the sole parameter determining the probability to 



find a LMXB in that GC or UCD (Sivakoff et al 



20071 ). The dependency of that probability is how- 



ever nevertheless almost linear to the encounter 
rate, while the dependency on th e metallicity is 
much weaker ( Jordan et al ] 120041 : ISivakoff et al 



20071 ). This may explain why the fraction of GCs 



with a LMXB is apparently already well explained 
if only the encounter rate in the GCs is considered 
(see Figure HJ despite the color-lu minosity rela- 
tion f or GCs in the Virgo cluster (cf. iMieske et al 



20061 ). It is thereby unlikely that the drastic dis- 



crepancy between the observed LMXB-frequency 
in UCDs and the theoretical prediction based on 
the encounter rate can be explained by an unac- 
counted metallicity effect, even though the color- 
luminosity dependency may be so mewhat more 



prono unced for UCDs than for GCs (jMieske et al 
I2OIOI) . 



The conclusion is that the large number of 
LMXBs in UCDs is best explained by a large num- 
ber dark remnants as a consequence of a top-heavy 
IMF in UCDs (and not as a consequence of differ- 
ent escape velocities or metallicities) . 

For an invariant IMF the theoretical LMXB 
frequency is highest at a luminosity of Ly « 
10^ Lq^, because in these systems the present- 
day stellar density has a maximum a nd close en- 



count ers are most frequent (Figure 4 in lDabringhausen et al 
I2OO8I ). If the very dense star formation conditions 
are responsible for a top-heavy IMF then, on first 
sight, the smallest IMF index is expected in 
systems with Ly ~ 10^ ^0,y and not in the most 
luminous UCDs. However, in systems with a top- 
heavy IMF stellar feedback is strongly enhanced 
and rapid gas expulsion leads to an exp ansion 
of the UCDs (|Dabringhausen et al.l lioToh . The 
UCDs revirialise after a few dynamical time scales 
(<100 Myr) and undergo no further size evolu- 
tion. Thus, their present day stellar density is the 
dynamically relevant quantity for producing the 
LMXB population. 
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We now determine by what amount the dark 
remnant content in UCDs has to be increased to 
get the theoretical LMXB-lrequency into agree- 
ment with the observed values. For this, Equa- 
tion (pS]) with 7 = 1 and 7 = 0.8 is used. 
This equation has and mtr as parameters 
(Equation [T]). In this paper, mtr — 1 Mq and 
mtr = 5 M0, so that the influence of the in 
principle quite arbitrary choice of mtr is tested. 
Note that with mtr — 1 M0, Equation (jlj de- 
sc ribes the family of IMFs th at were considered 



Dabringhausen et al.l (l2009[). For either choice 



of mtr, the canonical IMF (jKroupal l200ll 120021 ) 

corresponds to as = 2.3 and a smaller value of as 
increases the fraction of high-mass stars and sub- 
sequent dark remnants. The a^ that can explain 
the discrepancy between P{Lv) (i.e. the func- 
tion describing the observed LMXB frequency in 
UCDs) and Th{Lv) (i.e. the theoretical expec- 
tation for the LMXB frequency in UCDs if their 
IMF was canonical) at a given Ly can be found 
by numerically solving Equation (j26p for a^ with 
a given value for mtr- 

The Lv dependence of as required to bring the 
model into agreement with the UCD data is plot- 
ted as the solid line in Figure [S] for mtr = 1 Mq 
and in Fig ([6]) for mtr — 5 M©. In either case, 
the most massive UCDs must have an extremely 
top-heavy IMF in order to explain their LMXB- 
excess. The higher mtr is, the more exotic the IMF 
of UCDs must be in order to explain the number 
of LMXBs in them. For a given value for mtr, 
it is on the other hand only of minor importance 
whether P{Lv) is proprotional to Fh or propor- 
tional to Fjj . 

For mtr — 1 Mq, the independent analy- 
sis in this paper leads the same top-heavy IMF 
as derived from the UCD mass-to-light ratios 
( Dabringhausen et al. 20091) . shown as the dotted 



line in Figure [S] Such a comparison is not mean- 
ingful for mtr = 5 Mq, since the shape thereby 
assumed for the IMF is different f rom t he IMF 
considered in [Dabringhausen et al. ( 20091 ). 

The most likely relations between a^ and 
\ogiQ{Lv) shown in Figures [S] and [HI are remark- 
ably close to a linear function. 



as = clogio(Lv) + d. 



(27) 



The best fitting parameters c and d have been 
determined from a least-squares fit to 48 sample 



values calculated from Equation (1261) . These are 
shown Table [TJ Probably the best model for the 
IMF in UCDs is calculated when /lmxb oc F^-^ 
and mtr = 1 M© are assumed. This is because ob- 
servations suggest a less-than-linea r dependency 



of /t. mxr on the encoun ter rate ([Jordan et al 
2004t iSivakoff et al.l I2OO7I) and assuming 



mtr > 

1 M© implies even more extreme deviations from 
the canonical IMF in high-mass UCDs while the 
IMF is remark ably invariant in open star clusters 
(lKroupall200lh . 

Figure m suggests that the value of /lmxb for 
the most luminous UCDs is of central importance 
for estimating the slope of P{Lv) (Equation [^S]) 
and thus for the as calculated from Equation (|26l) . 
This is because of the distance of these data points 
to the other data points, which is due to the fact 
that the corresponding Ly interval is large. In or- 
der to estimate an uncertainty to the dependency 
of as on Ly , we changed the value of /lmxb for the 
most luminous UCDs [Ly > 2x 10® L©) by 3 times 
its uncertainty. P{Ly) was then recalculated with 
this new value and used in Equation (|26|). The 
resulting limits on the dependency of as on Ly 
are indicated by the limits to the gray area in Fig- 
ures [5] and [51 Also the limits of the gray areas 
are parametrized with linear functions, which are 
listed in Table [H 

The uncertainty of the upper mass limit for 
stars, mmax*, has little effect on the results sum- 
marized in Table [H This is illustrated with 
Figure |7l where the dependency between Ly 
and as calculated fron i Equation (l26l) is shown 
for mmax* = 150 M © (|Weidner fc Kroupalliool 



Oev fc Clarkd I2005D a nd for 



(jCrowther et al.l |2010[ ). The two 



= 300 M© 
functions are 



almost identical. 



4. The supernova rate in Arp 220 

A top-heavy IMF in UCDs can theoretically 
be understood if UCDs formed as very massive 
star clusters that were internally heated by infra- 
red radiation that was trapped inside a molecular 
cloud massive enough to form a UCD-type star 
cluster (|Murrav or if UCDs formed from 

molecular clouds that were heated by highly ener- 
getic cosmic rays originating from numerous type- 
II supernovae surrou nding those molecular clouds 
( Papadopoulod 2OI0I : cf. Section [T|). Both sce- 
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Fig. 5. — The IMF in UCDs for — 1 Mq. Plotted is the high-mass IMF index, a^, as a function of the V- 
band luminosity of the UCDs, Ly for /lmxb oh Fh leading to P{Lv) oc F^ (left panel) and for /lmxb cx: F°-® 

leading to P{Lv) oc Fj^' (right panel). The solid line shows the most likely high- mass index required to 
increase the dark remnant content in UCDs in order to match to observed LMXB-frequency (derived from 
the solid line in Figure 14]) . The grey shaded area marks an estimate for the 3 a region. The horizontal long 
dashed line marks the canonical IMF with = 2.3. The dotted line shows t he independently calculate d 
high-mass IMF index obtained from the observed mass-to-light ratios of UCDs ( Dabringhausen et al. 2009() . 
Simple-to-usc fitting relations for the variation of with Ly can be found in Tab. ([1]). 



narios imply that UCDs are formed during star- 
bursts, either because of the link between the 
formation of the most mas sive star-clusters an d 
high star formation rates (jWeidner et al 
or because the cosmic-ray field would only then 
be intense enough for effective heating of the 
molecular clouds. Note that likely progenitors of 
UCDs have actually been observed in star-bursts 
(|Fellhauer fc Kroupal[2Q02l ). 

Ultra-luminous infra-red galaxies (ULIRGs) are 
believed to be galaxi es with star-bursting regions 



(jCondon et al.l Il99ll ). They are thus systems 
where UCDs are probably forming. If this no- 
tion is correct and the IMF in UCDs is top-heavy, 
the ULIRGs as a whole should have more mas- 
sive stars than expected for an invariant, canon- 
ical IMF. As a consequence, the rate of type II 
supernovae is expected to be higher. 

In the following, we test the hypothesis of a 
top-heavy IMF in ULIRGs. For this reason, we 
quantify how the type-II supernova rate (SNR) in 
a star burst is connected to the star formation 



in it. Based on this, theoretical predictions for 
the SNR of Arp 220, which i s one of the closest 
ULIRGs I Lonsdale et al ] |2006l) . are calculated and 
compared to observations of this stellar system. 

The type-II supernova rate (SNR) observed in 
a stellar system depends on its IMF as well as on 
its star formation history (SFH), i.e. how the star 
formation rate in the stellar system has changed 
with time, because these quantities determine the 
numbers and ages of stars in given mass intervals. 
If star formation begins at a time toj only stars 
above a time-dependent mass-limit miow can have 
completed their evolution at a time t > to- For 
stars ev olving into SNe, this mass c an be approx- 
imated ( Dabringhausen et al. 2010[ ) by 



M^ 



= 74.6 



Myr 



- 2.59 



-0.63 



(28) 



Note that no stars evolve to type-II supernovae 
(SNe), if t-to < 2.59 Myr. 

Now consider a time interval [t, t + At] and stars 
in a mass interval [m,m + Am], where m > miow 
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fLMXB proportional to r 1^1^x6 proportional to r° ' 




1 2 4 6 10 1 2 4 6 10 

Lv [IO^Lq] Lv [lO^m 



Fig. 6.— The IMF in UCDs for mtr = 5 Mq, otherwise as Figure [S] 



If the SFR was constant for all t > to, the num- 
ber of stars evolving into SNe in the given mass 
interval during the time At is equal to the num- 
ber of new stars that arc formed in the same mass 
interval. Thus, 



ASNR 



SFR 



m+Am 



_i , am) dm (29) 

yr M© yr i 

in this case, where ^(m) is assumed to be given by 
Equation ([T|) with the normalization defined by 
Equation ([S]). This normalization keeps the total 
mass of the stars which are formed per unit time 
constant. 

If At is small compared to the time scale on 
which miow changes, the number of all stars that 
evolve during At can be approximated as 

SNR SFR , , 

^(m) dm. (30) 



yr 1 M© yr 1 

Note that the SFR in Equations ^ and (pO|) 
should be considered an average value over a time- 
scale t — tQ. Variations of the SFR on much shorter 
time-scales are of no importance here. 

The SFR of a ultra-luminous infra-red galaxy 
(ULIRG) can be estimated as 

SFR ipiR 



M© yr- 



5.8 X 109 Lr. 



(31) 



where Lpm is the far infra-red (FIR) luminosity 
of the ULIRG (|Kennicuttlll998h . 



One of the nearest ULIRGs is A rp 220. Using 
L fir = 1.41 X 10^2 for Arp 220 (jSanders et al 



20031 ) ■ Equation dSJ) implies a SFR of « 240 M, 



for that galaxy. The SNe in Arp 220 have been 
observed in a central region with a diameter of 
w 1 kpc, from where a bout 40 per cent o f its FIR 
luminosity originates ( Soifer et al. 19991 ). Equa- 
tion dsn thus implies a SFR of ^ 100 M©yr-i 
if only this part of Arp 220 is considered. Note 
that this SFR is consistent with the SFR that 
has been suggested for a forming UCD if UCDs 
form on a timescale of a pproximately 1 Myr 
([Dabringhausen et al. I I2OO9I) . Also note that the 
observed SN in Arp 220 do not seem to distributed 
evenly over the central part of Arp 220, but to be 
concentrate d in two knots which have a radius « 
50 pc each ( Lonsdale et al. 20061 ) (i.e. the size of 
a UCD). This implies that indeed a major part of 
the star formation in the central part of Arp 220 
takes place within these two knots. This would 
imply projected star formation densities of a few 
10"^ M© yr"i pc"^ in the knots. 

SNRs calculated from Equation ([30|) for a con- 
stant SFR of 100 M© yr^^ are shown as functions 
of the high-mass slope of the IMF in Figure [51 The 
two curves correspond to different times at which 
the star burst was initialized, but the expected 
number of SN per year (i.e. the SNR) is low in 
any case. The number of SN that actually occur 
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m.|.=1 Wlp m*|.=5 M,; 




0.8 ^ ' ' ' — 1 l-J ' ' ' — 1 

1 2 4 6 10 1 2 4 6 10 

Lv [IO^Lq] Lv [IO^Lq] 

Fig. 7. — The high- mass IMF slope, a^, as a function of F-band luminosity for different upper mass limits of 
the IMF, m„iax* assuming Wtr = 1 Mq (left panel) or Wtr = 5 Mq (right panel). The solutions, a^^Ly), are 
calculated from Equation (I^B)) with 7 = 1, i.e under the assumption that /lmxb oc Fh. Assuming 7 — 0.8 
leads however qualitatively to the same results as assuming 7 = 1 (c.f. Figures [5] and |6]). The solid line 
corresponds in both panels to mmax* = 150 Mq and is thus identical with the solid line in Figure [5l and 
Figure El respectively. The dashed lines corresponds to Wmax* — 300 Mq . 



within one year can therefore differ substantially 
from the calculated SNR, as the frequency of SN 
over such a time span obeys low-number statis- 
tics. Thus, the probability for a certain number of 
SN to happen within one year is quantified by the 
Poisson distribution function. 

Now consider the case that the star-burst in 
Arp 220 already lasts for more than 40 Myr. This 
implies rniow = 8 Mq , so that the number of SNII 
per year is maximized for the given SFR. The ex- 
pectation value for the SNII-rate is then about one 
per year if the IMF was canonical (i.e. = 2.3), 
but about two per year for a top-heavy IMF with 
1 ^ Q;3 ^ 2, where the SN-rate is only a weak 
function of (cf. Figure [5]). Thus, the probabil- 
ity to actually observe fou r new SNII in a given 
year ( Lonsdale et al. 20061 ) is then about two per 
cent if the IMF is canonical, but about 12 per cent 
for 1 < ag < 2. 

A more elaborate discussion of the SNII rate in 
Arp 220 is obtained by taking into account that 
stars in a galaxy form in star-clusters of differ- 
ent masses, since mmax of the IMF depends on 
the mass of the star-cluster for low-mass star- 



clusters. This implies that the integrated galactic 
IMF (IGIMF) of all star-clusters in Arp 220 com- 
bined is not equal to the IMF in its star-clusters. 
This IGIMF is given by 



IGIMF 



(to) 



<(SFR) 



^(m < m,„ax(Mocl)) 



X Cccl(Mccl)rfMccl, 



(32) 



where to is the initial stellar mass, Mcci is the 
initial stellar mass of a star cluster, Mcci,min is 
the minimum mass of newly formed star-clusters, 
Afeci,max(SFR) is the SFR-dependent maximum 
mass of newly formed star-clusters, ^{m) is the 
IMF and_£g£i(M££]j is the star-cluster mass func- 



tion ( Weidner &: Kroupa 2005 : IWeidner et al.l 
201 ll ). The IGIMF can be parametrized by a 
multi-power law. 



6GiMF(m) = kkim 



(33) 



with 



ai = 1.3, 
a2 = 2.3, 



O.K — < 0.5, 

- Mq 

TO 

0-5 < < 1, 

- Mq 
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Table 1: The best fitting parameters lor linear fits to P and lor different models. The different cases 
(most likely case, upper limit, lower limit) listed here lor every model correspond to different values ol 
P lor the UCDs with the highest masses (cl. Sec 13. 3p . Probably the best model lor the IMF in UCDs 
is calculated when /lmxb oc T^'^ and mtr = 1 Mq are assui ned. This is because observations sugg est a 



less-than- linear dependency ol /lmxb on the encounter rate ([Jordan et al.l l2004t ISivakoff et al.l 120071 ) and 
assuming mtr > 1 Mq implies even more extreme deviation s from the ca nonical IMF in high-mass UCDs 
while the IMF is remarkably invariant in open star clusters (jKroupall200lh . The parameters describing the 
IMF according to this model are shown in bold lace in this table. 











logio(^') 


= alogio(Ly) + b 


— c log 


',wiLv)+d 


model 








a 


b 


c 


d 


/lmxb 


oc 




: 1 Mq, most likely case 


0.207 


1.249 


-1.337 


2.332 


/lmxb 


oc 


Th, mtr = 


: 1 Mq, upper limit 


0.615 


1.201 


-1.878 


2.375 


/lmxb 


oc 


Th, mtr = 


= 1 Mq, lower limit 


-0.202 


1.298 


-0.884 


2.396 


/lmxb 


oc 


mtr 


= 1 Mq, most likely case 


0.207 


1.249 


- 1.402 


2.337 


/lmxb 


oc 


mtr 


— 1 Mq, upper limit 


0.615 


1.201 


2.089 


2.391 


/lmxb 


oc 


mtr 


= 1 Mq, lower limit 


-0.202 


1.298 


0.861 


2.304 


/lmxb 


oc 


Fh, mtr = 


: 5 Mq, most likely case 


0.207 


1.249 


-2.415 


2.275 


/lmxb 


oc 


Fh, mtr = 


: 5 Mq, upper limit 


0.615 


1.201 


-3.169 


2.289 


/lmxb 


oc 


Th, mtr = 


: 5 Mq, lower limit 


-0.202 


1.298 


-1.679 


2.263 


/lmxb 


oc 


^ h ' '"'tl' 


= 5 Mq, most likely case 


0.207 


1.249 


-2.512 


2.277 


/lmxb 


oc 


i h J '"-tr 


= 5 Mq, upper limit 


0.615 


1.201 


-3.442 


2.290 


/lmxb 


oc 


pO.8 


= 5 Mq, lower limit 


-0.202 


1.298 


-1.594 


2.263 



oigimf S 



1 < 



m 



M 



© 



where the factors ki ensure that the IGIMF is 
continuous where the power changes and fc is a 
normalization constant. Cigimf(to) equals il 
m < 0.1 Mq or ?Ti > Wmax*, where mmax* is the 
maximum stellar mass. Thus, the IGIMF defined 
here is equal to the IMF defined by Equation ([1]) , 
except lor the high-mass slope and the upper mass 
limit. 

The case ol a canonic al IMF in all star-clusters , 
2.3, implies (|Weidner fc Kroupal 120051 ) 



I.e. as 

CHGIMF ^ 3. The expectation value lor the number 
ol SNII per year would then be < 0.2 per year. On 
the other hand, aiciMF ^ 2 is possible, il a vary- 
ing IMF that becomes more top-heavy with star- 



cluster mass is considered (Weidner et al. 2011 



iKroupa et al. 2011 ). This implies that the proba- 
bility to_a£tuall};_oteervB^^ new SNII in a given 
year (jLonsdale et al is essentially zero il the 

IMF is canonical in all star-clusters, but it can still 
be about 10 per cent il the IMF becomes top-heavy 
in massive UCD-type star-clusters. 

The remnants produced by SNII are neutron 
stars and black holes. The SNII-rates thereby 



are an indicator lor how many mergers ol such 
remnants can be detected by searching lor grav- 
itational waves. Comparing the SN-rate lor 
aiGiMF = 3 to the SN-rate for oigimf — 2 thus 
suggests that about an order ol magnitude more 
ol such events may be expected il the IMF in 
massive star-clusters is not canonical, but top- 
heavy. Thus, the hitherto predicted detection rate 
ol about 30 mergers ol dark remnants per year 
( Baneriee et al. 2010l) lor the upcoming adLIGO- 
experiment could be too low by an order ol mag- 
nitude, as an invariant IMF has been used lor this 
estimate. 

Further evidence lor a top-hea vy IMF in star- 



bursti ng galaxies is found by [Anderson et al 
(|201l[ ) in Arp 299. They study numbers ol differ- 
ent types ol supernovae in Arp 299 and conclude 
from the mass ol the appropriate progenitor stars 
that the IM F is probably to p -heav y in that sys- 
tem. Thus, [Anderson et al. (|201ll ) qualitatively 
come to the same conclusion lor Arp 299 as we 
did lor Arp 220, while their method is different. 
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Fig. 8.— The SN-rate in the center of Arp 220. 
The SN-rates are functions of the slope of the IMF 
above a stellar mass of 1 Mq, a^, but also depend 
on the length of the star burst (which is indicated 
above the corresponding curve). They do not ex- 
ceed « 1 SN yr"-'^ for the canonical IMF (whose 
high-mass slope is marked by the dotted vertical 
line), or « 2 SN yr^^ for a top-heavy IMF. 



5. Star formation densities and the IMF 

A top-heavy IMF in UCDs is in-line with dif- 
ferent studies concluding a top-h eavy IMF in 
high- r edshift star formin g galaxies (jvan Dokkum 
20081 : iLoewensteinl 120061 ). Contrary to this, a 
recent spectroscopic study of two low-redshift 



very massive elliptical galaxies suggests a hith- 
erto unseen large popu l ation of low-mass stars 
(|van Do kkum k Conrovl l2010f l. which has been 
predicted as a possible cons equence of cooling 
flows on massive ellipticals ( Kroupa fc Gilmorg 
19941 ). It is on the other hand unlikely that the 
majority of the UCDs formed in potential wells 
deeps enough to cause cooling flows. 

Also note that the current stellar densities sug- 
gest that the star- formation densities (SFDs), i.e. 
the SFR per volume, of UCDs were very differ- 
ent from the SFDs of elliptical galaxies. Consider 
for instance an exemplar present-day UCD with 
AI — 10^ Mq and rh = 10 pc and an exemplar 
present-day elliptical galaxy with M = 10^^ M© 
and Th = 10** pc. These values can be consid- 



ered representative for typical UCDs and mas- 
sive elliptical galaxies, respectively (cf. figure 4 in 
Misgeld fc Hi lkcr 2011). Star formation is thought 



to have proceeded quickly in UCDs and massive 
elliptical galaxies, so that the stellar population 
of the exemplar UCD may have formed within 
lO'' yr ( Dabringhausen et al. 20091 ) and the stel- 
lar population of the exemplar el hptical galaxy 



may have formed within 10 yr ([Thomas et al 



20051 ). This leads to a SFR of 1 Mgyr"! for 
the exemplar UCD and to a SFR of 10^ Mg yr"! 
for the exemplar elliptical galaxy. The SFD can 
be estimated by dividing the SFR by rj^, lead- 
ing to a SFD of 10"^ Mgyr-^pc"^ for the ex- 
emplar UCD and a SFD of lO'^ Mq yr^^ pc^^ for 
the ex emplar elliptical galax y . How ever, accord- 
ing to Dabringhausen et al.l ( 2010t ) UCDs must 
have been even more compact when they formed 
(rh ~ 1 pc), since the mass loss following star 
formation with a top-heavy IMF must have ex- 
panded them to their present-day radii. With 
the masses of UCDs being 10^ Mq < Af < 
10® Mq, their SFRs ranged from 0.1 Mgyr^^ to 
10 Mq yr~^ if they formed within 10 Myr. An ini- 
tial rh of 1 pc thereby implies SFDs ranging from 
0.1 Mgyr-ipc-^ to 10 Mgyr-^pc-^. Thus, the 
SFDs of UCDs can easily be higher by six to ten 
orders of magnitude than the SFDs of massive el- 
liptical galaxies. Massive ellipticals can therefore 
not serve as a proxy for the stellar population in 
UCDs. 

It is therefore perhaps the SFD that deter- 
mines whether the IMF in some region of space 
becomes top-heavy, and not the overall SFR in 
a forming stellar system. This is actually consis- 
tent w ith models why the IMF ma y become top- 
heavy: [Dabringhausen et al. I (l2010l) argue that the 
central densities in forming UCDs were so high 
{p > 10^ Mqpc^^) that collisions and perhaps 
mergers between pre-stellar cores were important 
in them, in contrast to less massive stellar sys- 
tems. Likewise, if the heating of molecular clouds 
by cosmic rays is t he process by which the IMF be- 
comes top-heavy ( Papadopoulos 2010f ). it is again 
not the number, but the number density of the 
surrounding massive stars that makes heating of 
the molecular cloud effective. 
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6. Conclusion 

The dynamical mass-to-light ratios of ultra 
comp act dwarf galaxies (UCDs) are surprisingly- 



high (iHasegan et al 



20081: iMieske et al 



■ .2005 
20081). 



Dabringhausen et al 



This findin 



plained by iDabringhausen et al 



was ex- 
(120091) with an 



IMF that has more massive stars th an the canon- 
ical IMF deduced by iKroupa ( 200ll) from resolved 
stellar populations in the Milky Way. The high 
mass-to-light ratio of UCDs is then a consequence 
of a large population of dark remnants (i.e. neu- 
tron stars and black holes) in them. 

These dark remnants become visible as X-ray 
sources if they accrete matter from a low-mass 
companion star. The rate at which low-mass X- 
ray binaries (LMXBs) are formed in globular clus- 
ters and UCDs scales with the number density 
of dark remnants (see Section I3.2.ip . Data on 
the fraction of UCDs that harbour a bright X- 



ray source (jSivakoff et al.l 120071) can therefore be 
used to confirm the presence of a large population 
of dark remnants in UCDs by a method that does 
not rely on the fact that dark remnants only in- 
crease the mass of a UCD, but not its luminos- 
ity. It is shown in this paper that LMXBs in 
UCDs are indeed up to 10 times more frequent 
than expected for an invariant, canonical IMF. 
The overabundance of LMXBs is used to quan- 
tify the dependence of the high-mass IMF-slope, 
Q!3, on the luminosity of UCDs. This function is 
essentially equal to the dependence between the 
luminosity of the UCDs and their az suggested in 
IDabringhausen et al.l ()2009f ) based on the mass-to- 
light ratios of UCDs (see Section [33|) . Note that 
the Lv of present-day GCs and UCDs is just one 
of many properties of such systems. Dependencies 
of a3 on their initial mass, initial den sity and their 
metallitcity are therefore discussed in lMarks et al.l 
(|2012[ ). 

UCDs can be understood as the most mas- 
sive star-clusters which only form at extremely 
high galaxy-wide star fo rmation rates (SFRs) 
(jWeidner fc Kroupal liool ) . Alternatively, UCDs 
could form by the merger of gravitationally bound 
systems of star clu sters as they are observed i n in- 
teracting galaxies ([Fellhauer fc Kroupall2002l) . In 
either case, the formation of UCDs would be con- 
nected to star-bursts. Given that ultra-luminous 
infra-red galaxies (ULIRGs) are interpreted as 



galax ies with star-bursting regions ([Condon et al 
199ll ). they should show indications of a top-heavy 
IMF as a consequence. The nearest ULIRG is 
Arp 220. We show that the observed rate of 
type II supernovae in this ULIRG is indeed highly 
improbable if the IMF is invariant, but not if the 
IMF is top-heavy (see Section SJ. 

There are thus three mutually consistent ar- 
guments for a top-heavy IMF in UCDs or more 
generally star-bursting systems. Together with 
the evidence for the formation of UCDs being 
connected to star bursts, these arguments imply 
that the IMF becom es top-heavy in star-bursts (cf. 
Weidner et al.ll2011l ). This finding stands in con- 
trast to the prevalent notion that the IMF is in- 
variant (lKrouT)all2(30ll l2002t iBastian et allboiot 
Kroupa et al.l 120111 ) and thereby has important 
implications. For instance, estimates of the SFR 
of a galaxy based on observations that are sensi- 
tive only to high-mass stars and the assumption of 
an invariant IMF (like Equation [311) are too high 
if the IMF actually is top-heavy. Consequently, 
estimates for the time scale on which the pop- 
ulation of low-mass star in that galaxy is built 
up until the gas of the galaxy is depleted become 
too short. Also the chemical evolution of galaxies 
is different if the IMF in them can become top- 
heavy, since the nuclear reactions that occur in 
a star mainly depend on its mass. This has im- 
plications on their content of me tals and plane- 
tary systems ( Ghezzi et al. 2010t ). Furthermore, 
as more dark remnants are formed if the IMF is 
top-heavy, more dark-remnant mergers and thus 
gravitational- wave emitters should be detected in 
this case. Finally, the dynamical evolution of star 
clusters critically depends o n the shape of the IMF 
( Dabringhausen et al.ll2010l ). 
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